Carter-like constants of the motion in Newtonian gravity and electrodynamics 
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For a test body orbiting an axisymmetric body in Newtonian gravitational theory with multipole 
moments Qe, (and for a charge in a non-relativistic orbit about a charge distribution with the same 
multipole moments) we show that there exists, in addition to the energy and angular momentum 
component along the symmetry axis, a conserved quantity analogous to the Carter constant of Kerr 
spacetimes in general relativity, if the odd-£ moments vanish, and the even-^ moments satisfy Q21 = 
m(Q2/m) . Strangely, this is precisely the relation among mass moments enforced by the no-hair 
theorems of rotating black holes. By contrast, if Newtonian gravity is supplemented by a multipolar 
gravitomagnetic field, whose leading term represents frame-dragging (or if the electrostatic field 
is supplemented by a multipolar magnetic field), we are unable to find an analogous Carter-like 
constant. This further highlights the very special nature of the Kerr geometry of general relativity. 

PACS numbers: 04.70.Bw, 04.20.Jb 



I. INTRODUCTION 

The Kerr geometry of the rotating black hole in gen- 
eral relativity has long been notable for its remarkable 
properties (see, for example [J 0)- Apart from be- 
ing the unique, stationary, axisymmetric, asymptotically 
flat vacuum solution of Einstein's equations with a non- 
singular event horizon [3|, Hj , it manifests additional deep 
symmetries that are still not fully understood. One ex- 
pression of this is the statement that the Kerr black hole 
has "no hair"; its exterior spacetime geometry is com- 
pletely determined by two parameters, mass to and spin 
parameter a (the angular momentum per unit mass). For 
example, it is possible to characterize the exterior space- 
time of Kerr by a set of electric, or mass multipoles, Qi, 
and magnetic, or current multipoles, Jg, where £ > is 
an angular harmonic index. The near baldness of the 
geometry is expressed efficiently by the equation [|| 



}t + IJe = m(ia) 



(1) 



with Qq being the mass to, and J\ = S = ma being 
the angular momentum, in units where G = c = 1. As 
a consequence, for example Q2 = —ma 2 , and Q21. = 
m(Q 2 /m) e . 

Another remarkable property of the Kerr geometry is 
the existence of an additional constant of the motion for 
the geodesic of a test body, beyond the trivially con- 
served rest mass /i of the particle, the energy E, and 
the component of angular momentum L z along the black 
hole's symmetry axis. The conservation of E and L z arise 
simply from the stationarity and axisymmetry of the ge- 
ometry. The additional constant of the motion, called 
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the "Carter" constant is unexpected, since the ge- 
ometry is not spherically symmetric. It is a combination 
of the total orbital angular momentum L 2 , which itself 
is not conserved for generically oriented orbits, compo- 
nents of the particle's linear momentum (squared), and 
other variables dependent upon the angular momentum 
parameter a and the angle 9 between the particle's posi- 
tion and the symmetry axis 0, [§] • The existence of this 
added constant is intimately linked to the precise rela- 
tionship among the moments that characterize the de- 
viations from spherical symmetry, expressed by Eq. (fTj). 
Mathematically, it owes its constancy to the existence of 
a so-called "Killing tensor" £ Q/ 3, a generalization of the 
Killing vectors £ Q that generate the time and azimuthal 
symmetries that lead to the conservation of E and L z [9[ . 

A consequence of the existence of the Carter constant 
is that it makes the geodesic equation in Kerr completely 
integrable, that is, solvable in terms of quadratures. In 
practical terms, this may be a useful tool in analyses 
of the orbits of small black holes around massive black 
holes. These are the extreme mass-ratio inspirals (EMRI) 
that are an important potential source of gravitational 
waves for the proposed space interferometer LISA. Con- 
siderable work has gone into trying to determine how 
the Carter constant evolves with time because of gravi- 
tational radiation reaction, and how that plays into the 
evolution of inclined, inspiralling orbits around rotating 
black holes 0, [H El El • 

In a recent paper, Flanagan and Hinderer [3] stud- 
ied the evolution of the Carter constant using a post- 
Newtonian multipole approach for the effects of radiation 
reaction. In that paper, they pointed out the unusual fact 
that a Carter-like constant can exist even in Newtonian 
gravitational theory. For an axisymmetric body charac- 
terized by a mass to and quadrupole moment Q2, and 
to first order in Q2, there are three conserved quantities: 
energy per unit mass E = v 2 /2 — U, angular momentum 



2 



per unit mass along the symmetry axis e, h z 
where h = x x v, and a Carter- like constant 



C 



h 2 + ^ 
m 



(e-vf 



2Gm 



h ■ e 



(2) 



where n = x/r. (Since we are in Newtonian gravity, the 
conservation of rest mass p, which links coordinate time 
and proper time, is not relevant.) They went on to show 
that, for motion in a field with any individual non-zero 
harmonic Qi, with £ > 2, no such constant exists. 

However, this raises the question, in Newtonian grav- 
itation, and for that matter in Maxwell electrostatics, 
does a Carter constant exist for general axisymmetric 
distributions of mass (or charge)? In this paper we point 
out that the answer is yes, if the moments satisfy the 
constraints 



Q2e+i — , 

Qu = m(Q 2 /m) t 



(3a) 
(3b) 



Strangely, these are precisely the relations among the 
electric moments imposed by the no-hair theorem (TTJ) for 
the Kerr geometry. 

In this paper we prove this unusual property of Newto- 
nian gravity (Sec. [TTJ) . If one augments Newtonian grav- 
ity with a gravitomagnetic field (or, in electrodynamics, 
if one adds a multipolar magnetic field), we are not able 
to find an analogous Carter constant, unless the electric 
moments satisfy the conditions above, and all the gravit- 
omagnetic moments (or magnetic moments) vanish (Sec. 
IIIII) . In Sec. HVI we discuss the implications of this result. 



II. CARTER-LIKE CONSTANT OF THE 
MOTION IN NEWTONIAN GRAVITY 

The Newtonian equation of motion of a test body or- 
biting a massive body with a density distribution p(x) 
is 



a = d 2 x/dt 2 = VU. 



where 



U = G 



r*' 



\x — X' 



(4) 



(5) 



where G is the Newtonian gravitational constant and 
p{x) is the mass density. After suitable replacements 
involving charges, this is also the non-relativistic equa- 
tion of motion of a charge in the electrostatic field of a 
charge distribution. 

We Taylor expand l/\x — x'\ in powers of x' and define 
a set of symmetric, trace-free (STF) multipole moments 
/W, defined by 



I< L > := / p(x')x^d 3 x' : 



(6) 



where L denotes an ^-dimensional multi-index ijk . 



and x 



denotes a product of £ components of x' , 
made trace-free. For example x' 



(2) 



X X 



r' 5 lJ /3. See [15[ for discussion and definitions of multi- 
pole fields and STF tensors. We then obtain the equation 
of motion 



-G- 



1=2 



(-1)* 

a 



I {L) d lL - , (7) 



where r = \x\, and d lL denotes £ + 1 products of par- 
tial derivatives, with contraction performed over the £ 
repeated indices. The mass and center of mass of the 
body are chosen such that m = I = J pd 3 x, and 
P = = / px^^x. 

We assume that the central body is axisymmetric 
about an axis characterized by a unit vector e. It then 
follows that for I > 2 must be proportional to an 
STF product of £ components of e, in other words, 



7< L > := Q e e^ 



(8) 



The quantity Qi defines the mass or electric £-pole mo- 
ment, and can be shown to be given by 



pixypt^x , (9) 



where Pt{z) is a Legendre polynomial, and z := e ■ n. 
We make use of the well-known fact that 



d 



iL 



Q = {-\f +1 (2£+l)\l 



ML) 



(10) 



Then, substituting Eqs. © and (jTOJ) into the equation 
of motion ([7]), and exploiting the identity 



e m n {iL) = 



(2* + 1)1! 



<Pl +1 (z) - e'P^z)] , (11) 



where prime denotes d/dz, we obtain the equation of mo- 
tion 



ran \ - 



GQ e 



[nP' l+1 (z) - eP' t {z)} . (12) 



t=2 



Contracting this equation with v, using the fact that 

v = rn + rh , dPi/dt = (n ■ e)P[ , (13) 

where an overdot denotes d/dt, and using the recursion 
relation P' t+l — zP[ + (£ + l)Pe, it is straightforward to 
find a constant of the motion 

1=2 

This of course is the energy per unit mass, conserved 
because the potential is stationary. Defining the angular 
momentum per unit mass h := x x v, we find that 



dt 



(15) 



1=2 



3 



Thus the component h e of the angular momentum along 
the body's symmetry axis is conserved, a consequence of 
axisymmetry. 

But are there other non-trivial conserved quantities? 
Given that the only vectors in the problem are x, v, h 
and e, the only interesting possibility is the norm h 2 , 
which itself is conserved only for spherical symmetry or 
for an orbit confined to the equatorial plane. Calculating 
its rate of change, we find 



2 dt 



, dh 
h ■ — 

dt 



1=2 



1=2 



(16) 



Pulling out the 1 — 2 term, inserting the relation (TT3|) 
between n and v, it can be shown that the £ = 2 term 
can be written as 

3G ^ 2 (n • e) [(« • e) - f(n • e)] 



d 
dt 
d_ 
dt 



(n • e) 2 



Q2 /Gm 



QiGm 

r( n ■ e )( v ' e ) 

to r z 
1 



■("•e) --(v-e) 
to \ r 2 



(17) 



of the motion is then given by 

^ 1,2 , ^2/ ^2 o ™(<32/™) £+1 



£=0 

(20) 

This generalizes Eq. (2.8) of [Ij]. 

What kind of Newtonian source satisfies the conditions 
of Eq. ([3])? Its potential takes the form 

e=i 

Using the generating function for Legendre polynomials, 
J2^L t £ Pe(z) = (1 — 2tz + z 2 ) -1 / 2 , we can rewrite U in 
the form 



U(x) = 



Gm/2 Gm/2 



\x — ae\ \x + ae\ 



(22) 



1=2 



where a := {Qi/m) 1 / 2 . If Q2/m is positive, which cor- 
responds to the prolate case, then this is the potential 
of two point sources, each of mass m/2, separated by a 
distance 2a along the e axis. In the oblate case, a is 
imaginary, so the potential, while real, does not have an 
interpretation in terms of simple point masses. 



III. NEWTONIAN GRAVITY PLUS 
GRAVITOMAGNETISM 



where we used the equation of motion (|12p to eliminate 
Gmn/r 2 , pulled out other time derivatives, and used the 
recursion relation P[ = zP' t+1 — (£+ \)Pt+\. Combining 
this with the remaining terms in the sum in Eq. p6[) and 
again using Eq. (TIB")) , as well as further identities satisfied 
by the Legendre polynomials, we obtain finally, 



d 
~dt 



1 1,2 i 1 <?2 , -, 2 Qz i n GQi 



1=2 



e=o 



■$>-l)^fe--^-2 

1=2 v 



(18) 



No amount of further manipulation can convert the ex- 
pression on the right hand side into a total time deriva- 
tive. Since the Legendre polynomials are linearly inde- 
pendent, then for arbitrary orbits, there will be a con- 
stant of the motion of the chosen form if and only if 



= — Q e - 2 , forall£> 2. 



(19) 



Since Qi = 0, this implies Eq. ([3]). The Carter constant 



We now add gravitomagnetism to Newtonian theory. 
This is a model of gravity that encapsulates some, though 
not all, of the features of linearized general relativity 
(see [la], for example). In electromagnetism, it is the 
non-relativistic limit of the Lorentz-Dirac equation. The 
equation of motion in this case takes the form 



1 



a = VU + -v x (V x A„) , 

c 

where A g is the gravitomagnetic potential, given by 
AG fp^s. 



(23) 



C J \x — X' 

oo 



where J^ pL ^> are STF current multipole moments, given, 
for a stationary axisymmetric body by J(.j r ie^ pL '] here 
Ji = S is the ordinary angular momentum (or the mag- 
netic moment) of the body. 

This generates a term to be added to the equation of 
motion {7} given by 

1=1 



4 



This gravitomagnetic term does not affect the energy, 
but it does generate a new conserved angular momentum 
given by 



T h AG X- 



Ml - (e ■ n) 2 ] 

(e + iy 



(26) 



However, in this case, we find no Carter-like constant 
in general. Our attempt at construction included using 
a linear combination of h 2 , L 2 and L z S. We find that 
a Carter-like constant exists if conditions are met, 
along with Jg = for all t, which is our original purely 
Newtonian case. If S ^ 0, and if we work only to first or- 
der in multipole moments, then we find a Carter constant 
given by 

C = h 2 + ^(ve) 2 -2^(n-e) 2 ~— he, (27) 
m r cr 

which is consistent with Eq. (2.23) of [l4|. The absence 
of a constant in the general gravitomagnetic case high- 
lights the very special nature of geodesic motion in the 
Kerr geometry. 



the property that it is singular on a line joining the two 
masses; the singularity represents a strut required to hold 
the two masses apart. Does the Bach-Weyl solution pos- 
sess a Carter-like constant, or equivalently, does it admit 
an appropriate non-trivial Killing tensor? 

In Newtonian gravity, does there exist a physically rea- 
sonable, oblate distribution of matter (with Q2 < 0) that 
satisfies Eq. ©? 

What happens when the equations are extended to 
higher post-Newtonian orders? The gravitomagnetic 
model considered above included only terms linear in 
S = ma. If corrections to the equations of motion of 
order S 2 are included, can a Carter constant be found at 
some PN order and at higher order in the moments [1 91] ? 
Considerations such as these may elucidate what it is 
about Kerr that makes a Carter constant possible, and 
may also give insights into how to parametrize spacetimes 
that are "not quite" Kerr, in order to see how observa- 
tions of EMRIs could be used to test general relativity in 
the strong- field regime [20j |. 
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IV. DISCUSSION 

In a purely Newtonian context, this result may not 
have much practical significance, given the contrived na- 
ture of the source that possesses a Carter constant. How- 
ever, it does raise a number of related questions. 

The first is purely mathematical: in general relativity, 
there is an analogue to the solution for the field of two 
point masses separated by a fixed distance. Known as 
the Bach-Weyl solution [UH^], it is an exact, static ax- 
isymmetric vacuum solution of Einstein's equation, with 
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